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Abstract 

Visual  and  quantitative  information  were  obtained 
on  heat  transfer  and  flow  in  a  branched-duct  test  section 
that  had  several  significant  features  of  an  internal  cool¬ 
ing  passage  of  a  turbine  blade.  The  objective  of  this 
study  was  to  generate  a  set  of  experimental  data  that 
could  be  used  to  validate  computer  codes  for  internal 
cooling  systems. 

Surface  heat  transfer  coefficients  and  entrance  flow 
conditions  were  measured  at  entrance  Reynolds  numbers 
of  45,000,  335,000,  and  726,000.  The  heat  transfer  data 
were  obtained  using  an  Inconel  heater  sheet  attached  to 
the  surface  and  coated  with  liquid  crystals.  Visual  and 
quantitative  flow  field  results  using  particle  image  velo- 
cimetry  were  also  obtained  for  a  plane  at  midchannel 
height  for  a  Reynolds  number  of  45,000.  The  flow  was 
seeded  with  polystyrene  particles  and  illuminated  by  a 
laser  light  sheet. 

Computational  results  were  determined  for  the 
same  configuration  and  at  matching  Reynolds  numbers; 
these  surface  heat  transfer  coefficients  and  flow  velocities 
were  computed  with  a  commercially  available  code.  The 
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experimental  and  computational  results  were  compeired. 
Although  some  general  trends  did  agree,  there  were 
inconsistencies  in  the  temperature  patterns  as  well  as  in 
the  numerical  results.  These  inconsistencies  strongly 
suggest  the  need  for  further  computational  studies  on 
complicated  geometries  such  as  the  one  studied. 

Nomenclature 

A  heated  area,  m^ 

h  heat  transfer  coefficient,  W/m^  K 

Q  heat  rate,  W 

r  recovery  factor 

T  temperature,  ®C 

U^  free-stream  velocity,  m/sec 

V  kinematic  viscosity,  m^/sec 

p  density,  g/cm^ 

wall  shear  stress,  N/m^ 

Subscripts: 

e  electrical  heat  input 

L  loss 

LC  liquid  crystal 
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R  recovery 
s  static 

t  total 

Introduction 

A  continuing  objective  in  gas  turbine  technology  is 
higher  engine  efficiency.  One  method  of  obtaining  higher 
efficiency  is  to  use  higher  engine  operating  temperatures. 
Such  higher  turbine  inlet  temperatures  often  require  the 
use  of  internal  cooling  in  the  blades  and  vanes  to  main¬ 
tain  surface  temperatures  th{it  the  material  can  with¬ 
stand  for  sustained  periods  of  time.  When  internal 
cooling  is  used,  it  is  important  that  it  be  used  in  the 
most  effective  manner.  Effective  design  of  internal  cool¬ 
ing  passages  requires  an  understanding  of  the  flow  and 
temperature  patterns  that  occur,  and  it  is  greatly 
enhanced  by  the  availability  of  dependable,  validated 
computer  codes. 

An  experimental  study  was  performed  to  obtain 
quantitative  visual  information  on  flow  and  heat  trans¬ 
fer  in  a  test  section  that  had  several  significant  features 
of  the  internal  passage  of  an  experimental  cooled  radial 
turbine  blade  described  in  Ref.  1.  The  objective  of  this 
study  was  to  generate  a  set  of  experimental  heat  transfer 
and  flow  field  data  that  could  be  used  to  validate  com¬ 
puter  codes  for  internal  cooling  systems.  In  recent  years, 
there  have  been  several  attempts  to  numerically  model 
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internal  flow  and  heat  transfer  with  computer  codes.  ’ 
Reference  2  points  out  that  it  is  the  combination  of  a 
complete  experimental  database  and  accurate,  efficient 
analytical  tools  that  will  ultimately  shorten  design  times 
and  optimize  future  turbines.  Reference  4  further  em¬ 
phasizes  the  importance  of  coordination  and  cooperation 
between  experimental  and  computational  resesu-chers  in 
aeropropulsion  systems.  Reference  5  reports  on  a  code 
prediction  for  a  rotating  square  duct  that  was  verified 
by  experimental  data.  Through  such  coordination,  the 
results  of  the  present  study  will  be  used  to  validate 
internal  cooling  codes  for  more  complicated  geometries. 

Entrance  velocity  profiles  and  free-stream  turbu¬ 
lence  intensity  were  measured  to  establish  upstream 
boundary  conditions  for  the  computations.  Surface  heat 
transfer  coefficients  were  measured  at  entrance  Reynolds 
numbers  of  45,000,  335,000,  and  726,000  (based  on  the 
entrance  hydraulic  diameter).  The  heat  transfer  data 
were  obtained  by  using  a  uniform  Inconel  heater  sheet 
attached  to  the  test  section  floor  and  coated  with  ther- 
mochromic  liquid  crystals.  Liquid  crystals,  by  virtue  of 
their  property  of  changing  color  with  temperature,  pro¬ 
vide  a  means  of  measuring  surface  temperature  and  visu¬ 
alizing  thermal  patterns  at  any  desired  location.  The 
local  heat  transfer  coefficients  were  calculated  by  meas¬ 
uring  the  electrical  input  to  the  heater  sheet  under  flow 


conditions  and  using  a  heat  balance.  Reference  6  de¬ 
scribes  the  method  in  more  detail. 

Quantitative  visual  flow  field  meeisurements  from 
particle  image  velocimetry  are  also  presented  for  a  plane 
at  midchannel  height  and  a  test  section  entrance 
Reynolds  number  of  45,000.  The  flow  field  was  seeded 
with  l-/xm  polystyrene  particles  and  illuminated  by  a 
laser  light  sheet.  This  flow  field  work  is  an  extension  of 
work  done  in  conjunction  with  the  University  of  Akron 
where  the  procedure  and  the  computerized  data  reduc¬ 
tion  method  were  developed.  Details  of  this  method  are 
given  in  Ref.  7.  The  present  report  includes  a  much 
larger  extension  of  the  flow  field  covered  in  Ref.  7  such 
that  it  can  be  compared  with  the  heat  transfer  results. 

Computational  results  for  the  same  configuration 
were  obtained  with  a  commercially  available  code.  This 
code  was  used  to  calculate  surface  heat  transfer  coeffi¬ 
cients  and  flow  field  velocities  at  Reynolds  numbers 
matching  some  of  the  experimental  results.  The  experi¬ 
mental  and  computational  results  were  compared. 
Although  some  general  trends  did  agree,  there  were 
inconsistencies  in  the  heat  transfer  patterns  and  magni¬ 
tudes.  These  inconsistencies  demonstrate  the  need  for 
further  development  of  computational  studies  on  compli¬ 
cated  geometries  such  as  the  one  investigated. 

Apparatus 

Figure  1  is  a  schematic  of  the  complete  test  facil¬ 
ity.  Room  temperature  air  was  drawn  through  a  tunnel 
by  a  vacuum  exhaust  system.  The  two-piece  test  section 
(Fig.  2)  was  made  of  clear  acrylic  plastic  and  had  a 
peissage  height  of  5.08  cm.  The  entrance  section  was 
15.24-cm  wide,  and  the  main  test  section  was  divided  by 
a  partition  into  two  passages  12.7-cm  wide.  One  passage 
was  open,  and  the  other  passage  contained  three  rows  of 
0.95-cm  diameter  pins  in  a  staggered  array.  These  pins 
were  spaced  2.54-cm  apart  lateredly  in  rows  also  spaced 
2.54-cm  apart.  This  configuration  captures  the  signifi¬ 
cant  features  of  the  cooled  radial  turbine  blade  shown  in 
Fig.  3:^  namely,  the  flow  split  created  by  a  partition, 
the  sharp  comer,  the  cylindrical  pins  for  turbulence  gen¬ 
eration,  and  leading-edge  stagnation  regions.  Entrance 
velocities  were  mezisured  by  the  differential  between  a 
wedl  static  pressure  tap  and  a  total  pressure  probe  at  the 
entrance.  Pressure,  temperature,  and  flow  data  were 
recorded  by  a  commercial  data-acquisition  system. 
Photographic  data  were  taken  with  one  or  more  video 
cameras. 

For  the  heat  transfer  tests,  a  0.025-mm  thick 
Inconel  sheet  with  a  known  uniform  electrical  resistemce 
was  attached  to  the  floor  of  both  the  entrance  and  main 
test  sections  with  a  double-coated  adhesive  material. 
Each  sheet  had  copper  bus  bars  spot  welded  to  it,  was 
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coated  with  black  lacquer,  then  spray  painted  with 
liquid  crystals.  A  cross  section  of  the  test  surface 
(Fig.  4)  shows  the  2.5-cm-thick  foam  insulation  that 
prevented  heat  loss  through  the  floor.  For  the  heat 
transfer  tests,  the  cylindrical  pins  were  made  of  copper 
in  order  to  simulate  an  actual  blade  where  the  surface 
heat  is  transferred  to  the  pins.  The  pins  were  not 
attached  to  the  surface  but  rested  on  the  surface  simply 
by  gravity. 

The  test  section  for  the  flow  visualization  tests  was 
identical  in  configuration  to  the  heat  transfer  test  sec¬ 
tion  but  without  the  liquid-crystal  heater  sheets  and 
without  the  foam  insulation.  Another  difference  is  that 
the  pins  were  made  of  clear  acrylic  plastic  instead  of 
copper.  This  allowed  the  liiser  light  sheet  to  partially 
pass  through  the  pins  to  illuminate  the  flow  field. 

Procedure 

Liquid-Crystal  Calibration 

The  liquid-crystaJ  layer  was  calibrated  for  color- 
temperature  correspondence  by  a  thermocouple  attached 
to  each  heated  Inconel  sheet.  Diirk  blue  corresponded  to 
the  highest  temperature;  blue,  green,  yellow,  and 
reddish-brown  corresponded  to  decreasing  levels  of  tem¬ 
perature.  Yellow  occurred  over  the  narrowest  tempera¬ 
ture  band  of  0.06  ®C  and  was  therefore  the  color  used 
for  temperature  measurement.  The  yellow  color  location 
could  be  changed  by  varying  the  electrical  heat  input  to 
the  sheet  under  flow  conditions.  During  calibration, 
when  the  yellow  color  was  at  the  same  location  as  the 
thermocouple,  the  temperature  reading  was  37.8  “C. 
This  was,  therefore,  the  calibration  temperature. 

Experimental  Measurements 

Entrance  Velocity  Profile  and  Turbulence  Inten¬ 

sity.  Velocity  profiles  were  measured  near  the  upstream 
end  of  the  entrance  test  section  by  a  surface  static 
pressure  tap  and  a  traversable  total  pressure  probe.  As 
measurements  were  made,  the  probe  was  moved  in  small 
increments  from  the  floor  surface  to  a  location  just 
above  the  midchannel  height  position.  It  was  physically 
impractical  to  measure  the  profiles  at  higher  positions. 
Because  the  entrance  passage  was  symmetrical  about  the 
midchannel  height  position,  it  was  assumed  that  the 
velocity  profile  from  midheight  to  the  test  section  ceiling 
was  identical  to  the  measured  profile  from  midheight  to 
the  floor.  Free-stream  turbulence  intensity  was  meeisured 
by  hot  wire  anemometry.  AU  the  measurements  were 
made  at  entrance  Reynolds  numbers  of  approximately 
45,000,  335,000,  and  726,000. 

Heat  Transfer.  Heat  transfer  tests  were  performed 
at  the  same  entrance  Reynolds  numbers  used  for  the 


other  measurements:  45,000,  335,000,  and  726,000. 
These  corresponded  to  a  range  of  velocities  from  9.14  to 
155.5  m/sec.  The  operating  procedure  weis  initiated  by 
bringing  the  tunnel  to  the  velocity  necessary  to  produce 
the  desired  Reynolds  number.  Electric  power  was  then 
supplied  to  both  entrance  and  main  test  section  Inconel 
sheets  at  a  matching  heat  flux,  raising  the  temperature 
of  the  liquid  crystal  to  the  point  at  which  color  patterns 
appeared.  The  location  of  the  yellow  color  lines  was 
varied  over  the  surface  by  changing  the  power  input 
while  the  free-stream  velocity  remained  the  same.  Be¬ 
tween  data  recordings,  typically  about  15  min  was  given 
to  allow  thermal  equilibrium  to  occur.  Data  were  taken 
by  four  charged-couple  device  (CCD)  color  video  cam¬ 
eras  mounted  high  above  the  test  section.  One  camera 
viewed  the  entire  entrance  section,  another  camera 
viewed  the  entire  main  test  section,  and  the  final  two 
cameras  gave  closeup  views  of  the  cylindrical  pin  area. 
All  four  cameras  were  activated  simultaneously  to 
achieve  time  correspondence,  and  the  data  were  recorded 
on  Super  VHS  video  tape.  A  photograph  of  the  test  sec¬ 
tion  under  flow  conditions  in  shown  in  Fig.  5.  The  color- 
temperature  patterns  on  the  surface  are  clearly  visible. 

Flow  Visualization.  Flow  visualization  tests  were 
performed  at  only  one  entrance  Reynolds  number, 
45,000.  In  these  tests,  the  flow  was  seeded  with  l-/im 
polystyrene  spheres  that  followed  the  flow  field.  These 
particles  were  suspended  in  a  heated  water-ethyl  alcohol 
mixture  and  sprayed  into  the  flow  stream  with  a  spray 
nozzle.  By  the  time  the  mixture  reached  the  test  section, 
much  of  the  liquid  had  evaporated.  The  particles  were 
illuminated  by  a  7-W  eu-gon-ion  laser  beam  that  was 
focused,  through  a  series  of  mirrors  and  lenses,  to  a  sheet 
of  light  approximately  2  mm  thick.  The  result  was  a 
thin,  illuminated,  horizontal  slice  of  the  flow  field  at 
midchannel  height. 

A  single  black  and  white  CCD  video  camera  re¬ 
corded  the  flow,  showing  the  trajectories  of  the  particles 
as  streaklines  on  each  video  frame.  A  light  chopper  was 
placed  in  the  laser  beam  to  measure  the  velocity  of  the 
flow  field.  This  chopper,  which  consisted  of  a  rotating 
disk  with  holes  spaced  such  that  the  light  wets  periodi- 
ceiUy  blocked,  could  chop  light  at  a  maximum  frequency 
of  3000  Hz.  The  instantaneous  velocity  of  the  particles 
was  represented  by  the  distance  between  interruptions  in 
the  streaklines  on  any  one  video  frame  as  illustrated  in 
Fig.  6.  Similarly,  the  acceleration  of  the  particles  was 
represented  by  the  change  in  the  distance  between  inter¬ 
ruptions  in  the  streaks.  Because  of  the  high  magnifica¬ 
tion  required  to  see  the  particle  traces,  only  one  small 
area  could  be  viewed  at  a  time.  A  computer-aided  design 
(CAD)  software  package  was  used  to  combine  the  data 
from  approximately  300  individual  frames  so  that  all  of 
the  information  could  be  presented  in  one  composite 
picture. 
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Data  Reduction 


Heat  Transfer.  Surface  heat  transfer  coefficients 

were  calculated  as  follows: 

A(Tlc  -  Tr) 

The  local  heat  transfer  coefficient  h  was  calculated  at 
the  location  of  the  calibrated  color  band  (an  isotherm — 
which  in  this  case  also  represents  a  uniform  heat  transfer 
coefficient).  The  heat  energy  supplied  to  the  heater 
sheet  was  calculated  from  the  measured  voltage  across 
the  sheet  and  the  current  through  a  shunt  resistor  in  ser¬ 
ies  with  the  sheet.  The  heat  loss  Ql  was  the  calculated 
radiation  from  the  surface.  The  insulation  under  the  test 
cell  floor  prevented  back  side  heat  loss,  and  other  heat 
losses  were  considered  negligible.  The  area  A  was  the 
measured  area  of  the  Inconel  heater  sheet,  and  the  tem¬ 
perature  of  the  test  surface  T^c  was  the  calibrated 
liquid-crystal  temperature.  The  recovery  temperature 
Tr  of  the  entrance  free-stream  air  was  calculated  as 

Tr  =  T,  +  r(Tt  -  T3) 

where  T^  and  Tj  are  the  static  and  total  air  tempera¬ 
tures,  respectively,  and  r  is  the  recovery  factor,  defined 
as  the  cube  root  of  the  Prandtl  number. 

An  uncertainty  analysis  based  on  the  method  de¬ 
scribed  in  Refs.  8  and  9  wm  performed  on  the  calcula¬ 
tion  of  experimental  heat  transfer  coefficients.  With  this 
method,  the  maximum  expected  uncertainty  values  in 
the  calculation  of  heat  transfer  coefficients  were  14.0, 
3.7,  and  3.8  percent  for  Reynolds  numbers  of  45,000, 
335,000,  and  726,000,  respectively.  The  relatively  high 
percentage  heat  transfer  uncertainty  at  the  lowest 
Reynolds  number  resulted  from  a  high  uncertainty  in 
flow  velocity  because  of  the  lower  accuracy  of  the 
pressure-measuring  devices  at  very  small  pressure  differ- 
entiab.  At  all  Reynolds  numbers,  errors  were  abo  possi¬ 
ble  in  reading  the  exact  location  of  the  botherms. 

The  heat  transfer  video  data  were  reduced  with  a 
personal  computer  equipped  with  a  commercial  video 
frame  grabbing  board.  A  computer  program  was  written 
to  manually  digitize  (trace)  the  yellow  color  bands  from 
a  video  image  grabbed  from  the  video  tape,  and  the 
digitbed  image  was  stored  in  a  computer.  The  digitized 
locations  along  with  the  measurements  of  surface  heat 
flux,  air  temperature,  and  free-stream  velocity  allowed 
the  calculation  of  heat  transfer  coefficients. 


Flow  Vbualization.  The  method  for  reducing  the 
flow  vbualization  data  b  described  in  detail  in  Ref.  7.  In 
the  present  study,  the  region  of  the  test  section  mapped 
with  thb  technique  was  greatly  increased  to  include  an 
area  from  8.9-cm  upstream  of  the  front  row  of  pins  to 
3.8-cm  downstream  of  the  last  row  of  pins. 

The  uncertainty  involved  in  determining  velocity 
is  largely  dependent  on  the  accuracy  of  the  operator  dur¬ 
ing  the  digitizing  process.  Reference  7  reports  that  maxi¬ 
mum  expected  errors  may  vary  from  4  to  8  percent. 

Computational  Solutions 

Computational  solutions  of  the  flow  field  and  heat 
transfer  coefficients  along  the  floor  of  the  test  section 
were  obtained  with  the  commercially  available  code, 
STAR-CD.^®’^^  This  code  solves  the  fully  three- 
dimensional,  turbulent  compressible  Navier-Stokes  equa¬ 
tions  on  an  unstructured  mesh  with  local  mesh  refine¬ 
ment  capabilities.  The  governing  equations  of  mass, 
momentum,  energy,  turbulent  kinetic  energy,  and  dbsi- 
pation  are  dbcretbed  by  the  finite  volume  method, 
and  a  semi-implicit  formulation  b  utUbed  (SIMPLE 
scheme). 

The  computational  mesh  for  this  analysb  b  shown 
in  Fig.  7.  Because  computer  resources  were  limited,  the 
computational  grid  was  limited  to  about  250,000  ceUs, 
with  finer  mesh  generated  in  the  pin  region.  Boundary 
conditions  were  set  at  the  entrance  plane  and  at  the  exit 
of  the  computational  domain  to  match  those  of  the 
experimental  values  at  Reynolds  numbers  of  45,000, 
335,000,  and  726,000.  Because  of  difficulties  encoun¬ 
tered  trying  to  input  the  measured  velocity  profile  into 
the  code,  a  constant  velocity  profile  was  used  at  the 
entrance.  This  profile  was  the  free-stream  velocity  meas¬ 
ured  at  the  entrance  plane  of  the  test  section.  In  addi¬ 
tion,  the  average  experimentally  measured  turbulence 
level  of  2.3  percent  was  specified.  A  constant  pressure 
was  assigned  at  the  exit,  and  the  walls  were  specified  to 
be  at  a  constant  temperature. 

Steady-state  calculations  with  the  upwind  differ¬ 
ence  scheme  were  chosen  for  this  analysb.  The  k-e  tur¬ 
bulence  model  was  used,  and  heat  transfer  coefficients 
were  obtained  with  the  Reynolds  analogy,  which  relates 
the  convective  heat  transfer  coefficient  to  the  shear 
stress  of  the  flow  at  the  wall.  A  Unix-based,  Silicon 
Graphics’  Irb  workstation  with  four  central  processing 
units  (cpu’s)  in  parallel  was  used  to  run  the  analysb  and 
to  pre-  and  post-process  the  data.  The  cpu  time  for 
obtaining  a  convergent  solution  ranged  from  29  to  55  hr. 
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The  converged  solution  of  the  flow  provided  veloci¬ 
ties,  pressures,  temperatures,  density  variations,  and 
various  other  parameters  throughout  the  channel,  as 
well  as  heat  transfer  coefficients  on  the  surface. 

Results  and  Discussion 
Experimental  Results 

Entrance  Velocity  Profiles.  Velocity  profiles  at  the 
entriince  are  shown  in  Fig.  8  for  the  three  entrance 
Reynolds  numbers  of  45,000,  335,000,  and  726,000. 
Fully  developed  flow  had  not  occurred,  and  as  expected, 
there  was  a  slight  decrease  in  boundary-layer  thickness 
as  velocity  was  increased.  The  boundary-layer  thickness, 
defined  as  the  point  where  the  velocity  is  99  percent  of 
the  free-stream  value,  was  1.55  cm  at  a  Reynolds  num¬ 
ber  of  45,000,  1.40  cm  at  335,000,  and  0.94  cm  at 
726,000. 

The  velocity  data  were  also  plotted  logarithmically 
in  terms  of  dimensionless  coordinates  (Fig.  9).  This  plot 
was  made  to  show  how  the  boundary  layer  compared 
with  turbulent  boundary  layers  normally  encountered  in 
turbine  flow.  Included  in  the  figure  is  the  calculated 
theoretical  line  for  a  fuUy  turbulent  boundary  layer.  The 
equations  used  were  taken  from  Ref.  13.  The  experimen¬ 
tal  data  match  the  theoretical  line  fairly  well. 

Free-Stream  Turbulence  Intensity.  Free-stream 
turbulence  intensity  weis  measured  with  a  hot  wire  ane¬ 
mometer  at  midchannel  height.  At  a  Reynolds  number 
of  45,000  the  turbulence  intensity  was  mesisured  as 

2.8  percent,  at  a  Reynolds  number  of  335,000  it  was 
1.2  percent,  and  at  a  Reynolds  number  of  726,000  it  was 

2.9  percent.  The  average  intensity  was  2.3  percent.  It 
has  not  been  determined  why  the  turbulence  intensity 
wM  so  much  lower  at  the  medium  value  of  Reynolds 
number. 

Heat  Transfer  Measurements.  Quantitative  heat 
transfer  coefficients  on  the  floor  of  the  test  section  were 
mapped  for  three  different  entrance  Reynolds  numbers. 
A  typiced  example  of  the  results  is  shown  in  Fig.  10. 
Parts  (a)  to  (c)  are  actual  video  screen  photographs  of 
the  surface  at  the  lowest  Reynolds  number  of  45,000  and 
at  three  of  the  various  electrical  energy  inputs.  The 
figure  includes  closeup  views  of  the  color-temperature 
patterns  around  the  pins.  Figure  11(a)  shows  the  corre¬ 
sponding  map  of  heat  transfer  coefficients  calculated 
from  these  isotherms,  and  Fig.  11(b)  is  a  closeup  of  the 
area  around  the  pins.  Similar  heat  transfer  maps  are 
shown  in  Figs.  12  and  13  for  Reynolds  numbers  of 
335,000  and  726,000,  respectively. 

Generally,  the  heat  transfer  patterns  were  similar 
for  all  Reynolds  numbers.  The  heat  transfer  in  the 


entrance  region  decreased  going  downstream  as  shown  in 
Figs.  11(a),  12(a),  and  13(a).  Areas  of  low  heat  transfer 
are  seen  in  the  front  left  (relative  to  the  flow  stream) 
corner  of  the  main  test  section  where  the  flow  is  recircu¬ 
lating.  The  lowest  values  of  heat  transfer  occurred  in 
this  region  at  all  three  flow  conditions  as  shown  in  these 
same  figures.  In  the  unobstructed  passage  (no  pins)  there 
is  a  similarity  in  the  pattern  of  the  heat  transfer  at  all 
three  flow  conditions  with  a  high  value  heat  transfer 
coefficient  coming  off  the  ptirtition  leading  edge  (proba¬ 
bly  produced  by  the  accelerating  flow  and  high  turbu¬ 
lence  created  in  this  area,  as  will  be  discussed  later).  The 
heat  transfer  coefficients  decreased  as  flow  progressed 
downstream.  Close  to  the  partition,  the  coefficients  were 
even  smaUer,  probably  because  of  a  flow  separation  in 
this  area.  From  the  closeup  views  of  the  pin  region 
(Figs.  11(b),  12(b),  and  13(b)),  the  cooling  effectiveness 
of  the  pins  is  evident  by  the  color  patterns  immediately 
downstream  of  the  pins  and  the  corresponding  high  val¬ 
ues  of  heat  transfer  coefficients.  Further  downstream  of 
the  pins,  the  heat  transfer  coefficients  decreased.  At 
some  locations,  the  experimental  contours  were  broad 
bands  instead  of  thin  lines,  indicating  very  gradual  sur¬ 
face  temperature  gradients. 

Flow  Field  Measurements.  Flow  field  measure¬ 
ments  were  made  at  an  entrance  velocity  of  9.14  m/sec, 
which  corresponded  to  the  lowest  Reynolds  number  of 
45,000.  Attempts  were  made  to  obtain  results  at  the 
higher  Reynolds  numbers,  but  these  attempts  were  un¬ 
successful  with  the  equipment  that  was  available.  As  the 
velocity  of  the  flow  field  weis  increased,  it  was  more  dif¬ 
ficult  to  detect  individual  traces  on  the  video  screen. 
The  maximum  measurable  particle  velocity  was  further 
limited  by  the  maximum  frequency  at  which  the  light 
chopper  could  operate,  which  was  3000  Hz.  As  previ¬ 
ously  shown  (Fig.  6),  the  velocity  was  determined  by  the 
distance  between  interruptions  in  the  streaks  shown  in 
any  one  video  frame.  At  higher  velocities,  streaks  would 
appear  in  a  single  frame  with  either  one  interruption  or 
no  interruptions  at  all,  making  it  impossible  to  deter¬ 
mine  the  velocities  or  accelerations. 

Figure  14(a)  is  a  video  picture  of  the  particles  as 
they  flowed  between  a  pin  and  the  partition  leading 
edge.  Figure  14(b)  shows  the  data  after  being  processed 
into  velocity  vectors.  The  numbers  on  the  diagram  iden¬ 
tify  each  particular  trajectory,  and  the  magnitude  of  the 
velocity  is  indicated  by  the  length  of  the  arrows. 

Figure  15  is  a  composite  view  of  about  300  video 
frames,  giving  an  overall  view  of  the  flow  field  in  the 
area  investigated.  Velocity  vector  diagrams  similar  to 
the  one  shown  in  Fig.  14(b)  can  be  generated  anywhere 
within  this  area  from  the  data  collected  in  this  study. 
Figure  15  reveals  a  few  features  that  help  explain  some 
of  the  heat  transfer  results.  A  turbulent  area  appears 
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immediately  downstream  of  the  pins,  which  corresponds 
to  the  high  heat  transfer  area  seen  in  the  maps.  In  the 
unobstructed  passage  where  the  flow  comes  off  the  parti¬ 
tion,  there  seems  to  be  a  region  of  turbulent  flow  bor¬ 
dered  by  an  area  of  smooth  accelerated  flow.  This 
corresponds  to  the  high  values  of  heat  transfer  seen  on 
the  heat  transfer  maps.  In  the  front  left  comer  of  the 
main  test  section,  the  flow  mapping  appears  to  indicate 
a  highly  turbulent,  high  velocity  area;  but  actual  obser¬ 
vation  of  the  flow  showed  this  area  to  have  very  slow, 
three-dimensional  flow,  making  it  diSlcult  to  properly 
visualize  the  flow  in  a  thin  slice  of  light.  The  seeded  air 
seemed  to  just  circle  slowly  around  in  a  lazy  manner, 
which  would  lead  to  the  low  heat  transfer  and  high  sur¬ 
face  temperatures  indicated  in  the  heat  transfer  tests. 

Computational  Results 

Heat  Transfer  Computations.  Computational  solu¬ 
tions  for  the  surface  heat  transfer  coefficients  were 
obtained  for  all  three  Reynolds  numbers.  The  results  are 
shown  in  Fig.  16  for  a  Reynolds  number  of  335,000  for 
comparison  with  the  experimental  results.  Figure  16(a) 
shows  a  decreasing  value  of  heat  transfer  coefficient  in 
the  entrance  section.  Low  values  were  seen  in  the  recir¬ 
culation  region  as  expected.  In  the  area  downstream  of 
the  pins,  the  isotherms  occurred  in  a  generally  longitudi¬ 
nal  direction.  Figure  16(b)  shows  a  closeup  view  of  the 
area  around  the  pins.  Low  values  were  calculated  in  the 
wake  of  the  pins. 

Flow  Field  Computations.  Computational  solutions 
for  velocities  at  midchannel  height  of  the  flow  field  are 
shown  in  Fig.  17.  The  color  chart  designates  the  magni¬ 
tude  of  the  velocities.  These  solutions  were  obtained  for 
a  Reynolds  number  of  45,000,  which  corresponds  to  an 
entrance  velocity  of  9.14  m/sec.  As  the  flow  approached 
the  pins  and  the  partition  area,  it  accelerated  into  the 
unobstructed  peissage,  whereas  the  flow  going  through 
the  pins  decelerated.  The  recirculation  region  shows  very 
low  velocities  in  the  range  of  2  to  4  m/sec.  The  flow  in 
the  pin  area  does  not  seem  to  be  as  turbulent  as  ex¬ 
pected.  Evidently,  the  grid  around  the  pins  weis  not  fine 
enough  to  show  the  complex  mixing  that  occurs  in  the 
wake  regions. 

Computational-Elxperimental  Comparisons 

Heat  Transfer  Comparisons.  There  is  agreement  in 
some  of  the  general  trends  between  the  experimental  and 
computational  results  (for  example,  the  decreasing  heat 
transfer  in  the  entrance  region),  but  there  is  disagree¬ 
ment  in  the  numerical  values  around  and  downstream  of 
the  pins.  Because  the  Reynolds  analogy  was  used  in  the 
analysis,  this  calculation  would  not  be  expected  to  be 
valid  around  the  pins.  In  addition,  the  computational 
map,  which  shows  longitudinal,  fingerlike  contours 


downstream  of  the  pins,  differs  from  the  exi)erimental 
maps,  which  showed  almost  constant  heat  transfer  lines 
across  the  cheinnel.  Figure  16  shows  high  heat  transfer 
coefficient  values  in  the  high  acceleration  region  between 
the  partition  emd  the  recirculation  area,  whereas  Figs.  11 
to  14  show  high  heat  transfer  coefficient  values  in  the 
high  turbulent  flow  coming  off  the  partition.  The  use  of 
a  finer  mesh  grid  near  the  floor,  where  the  heat  transfer 
is  calculated,  might  have  produced  more  compatible 
results.  In  the  recirculation  region,  the  calculation  shows 
a  low  heat  transfer  coefficient,  which  is  confirmed  by  the 
experiment. 

Flow  Field  Comparisons.  Because  experimentad 
flow  field  data  were  obtained  at  only  one  entrance  flow 
velocity,  9.14  m/sec,  a  comparison  between  computa¬ 
tional  aind  experimental  results  could  be  made  only  at 
the  one  corresponding  Reynolds  number  of  45,000  and 
only  at  one  small  area  at  a  time.  The  area  chosen  was 
the  region  near  the  leading  edge  of  the  partition  and  the 
nearest  pin.  A  close-up  view  of  the  flow  in  this  area  is 
shown  in  Fig.  18.  (This  view  is  at  midchannel  height.) 
In  the  area  upstream  of  the  pins,  there  b  good  agree¬ 
ment  with  experimental  velocities  at  about  7.5  m/sec. 
However,  in  the  area  between  the  partition  and  the  pin, 
the  code  showed  the  flow  decelerating  to  5.5  m/sec, 
which  b  lower  than  experimental  values  seen  in 
Fig.  14(b).  The  calculated  values  for  the  flow  field  show 
the  general  trends,  but  failed  to  predict  the  leveb  and 
details  around  the  pins.  The  grid  that  was  used  to  com¬ 
pute  the  flow  field  around  the  pins  consbted  of  16  celb 
around  the  pins  and  20  cells  deep — evenly  spaced  from 
top  to  the  bottom  of  the  channel.  The  number  of  the 
grid  points  was  chosen  arbitrarily  to  maintain  a  reason¬ 
able  size  grid  and  to  obtadn  a  solution  in  a  timely 
mzinner.  Using  a  denser  grid  or  assuming  symmetry 
about  the  midchannel  height  and  analyzing  only  the 
lower  half  with  the  same  number  of  grid  points  might 
have  generated  a  better  flow  field  calculation.  Such 
refinements  will  be  attempted  at  a  later  date. 

Concluding  Remarks 

The  experimental  investigation  revealed  the  tem¬ 
perature  patterns  that  occur  throughout  a  test  section 
with  some  significant  features  of  an  internal  passage  of 
a  cooled  turbine  blade.  Areas  of  high  and  low  heat 
transfer  were  identified,  and  quantitative  heat  transfer 
data  were  obtained  over  a  range  of  flow  conditions.  The 
cooling  effectiveness  of  the  cylindrical  pins  was  demon¬ 
strated.  Areas  of  recirculating  and  sepeirating  flows  and 
their  effect  on  surfewe  temperature  patterns  were  evi¬ 
dent.  The  experimental  flow  vbualbation  data  generally 
supported  the  heat  transfer  data  in  terms  of  explaining 
some  of  the  temperature  patterns.  Two  benchmark  sets 
of  data  were  produced,  heat  transfer  and  flow  field. 
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which  are  now  available  for  interned  cooling  code 
validation. 

Computational  solutions  for  the  flow  field  and 
surface  heat  transfer  were  obtained  with  a  commercially 
available  code,  STAR-CD.  These  solutions  showed  gen¬ 
eral  trends  of  the  flow  field  in  the  channel.  The  mag¬ 
nitude  of  the  velocity  field  around  the  pins  was  lower 
than  the  experimental  vedues.  The  code-predicted  heat 
transfer  results  agreed  with  some  of  the  general  trends  of 
the  experimental  data,  but  there  were  differences  in  the 
numerical  values  of  the  surface  heat  transfer  coefficients 
on  most  of  the  channel  floor.  In  the  pin  area,  the  calcu¬ 
lated  heat  transfer  coefficients  downstream  of  the  pins 
showed  trends  opposite  to  those  of  the  experiment.  The 
grid  size  of  250,000  cells  was  evidently  insufficient  to 
resolve  the  flow  field  in  this  area  as  well  as  the  heat 
transfer. 

STAR-CD,  which  is  typical  of  currently  available, 
general  purpose,  computational  fluid  dynamics  (CFD) 
codes,  showed  the  limitations  that  exist  at  the  present 
time  in  predicting  the  flow  field  and  heat  transfer  in 
complex  geometries.  The  difficulties  encountered  in 
attempting  to  match  computational  simulation  and  ex¬ 
perimental  data  demonstrate  the  need  for  further  devel¬ 
opment  of  CFD  codes. 
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Figure  1 . — Branched  duct  tunnel.  Note:  Only  one  camera  used  for  flow-visueilization  tests. 
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Figure  3. — Cooled  radial  turbine  blade  (Ref.  1). 
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Figure  5— Test  section  underflow  conditions;  Reynoids  number,  45,000. 
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Figure  6. — Instantaneous  velocity  from  1 /30-sec  video  frame 
using  3000-Hz  light  chopper. 


Figure  7. — Computational  mesh. 


13 


Velocity,  u*  Distance  from  floor  (cm) 


5 


4 


o  Re  =  45,000 

X  Re  =  335,000 

a  Re  =  726,000 

(5  Boundary  Layer 

Thickness 


3  - 


o 

o 

-O- 


X 

X 


□ 

□ 


"X" 

X 


X 

X- 

X 


Mid-channel 


o 

o 

o 

o 

o 


o 

o 

o 

o 

o 

o 

o 

o 

0 

o 

jn  X. 


V' 


X 

Lx 


8 

□ 

□ 

□ 

□ 

□ 

□ 

□ 


q_ 

a 

□ 


m 


20  40  60  80  100  120  140  160  1 80  200 

Velocity  (m/sec) 


Figure  8.— Entrance  velocity  profiles. 
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Figure  9. — Dimensionless  boundary-layer  profile. 
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(a)  Low-heat  input;  main  section. 

Figure  10. — Surface  coior-temperature  patterns;  Reynolds  number  45,000. 
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(b)  Medium  heat  input. 
Figure  1 0. — Continued. 
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(c)  High  heat  input. 

Figure  1 0. — Conciuded. 

I 


19 


3AH. 


Page  intentionally  left  blank 


Separation  region 


(a)  Entire  test  surface. 


(b)  Close-up  of  pin  area. 


Figure  1 1  .—Experimental  heat  transfer  coefficients;  Reynolds  number,  45,000.  (Coefficients  given  in  W/m^.K.) 
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Separation  region 


(b)  Close-up  of  pin  area. 


Figure  1 2.— Experimental  heat  transfer  coefficients;  Reynolds  number,  335,000.  (Coefficients  given  in  W/m^.K.) 


Separation  region 


(a)  Entire  test  surface. 


(b)  Close-up  of  pin  area. 


Figure  13. — Experimental  heat  transfer  coefficients;  Reynolds  number,  726,000.  (Coefficients  given  in  W/m2.K.) 
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(a)  Video  picture. 


Figure  14. — Particle  trajectories;  Reynolds  number,  45,000. 
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Figure  1 5. — Composite  of  digitized  video  data;  Reynolds  number,  45,000.  (300  video  frames  combined.) 


25 


Recirculation  region 


26 


p  Recirculation  region 

J _ 

\ 

\ 

\ 

\ 


c 


1  O- 

0  . 

.  .....  V  W  ...  .  ,  j 

•  1 

1 

o 

~  ■  1 

:  o 

o 

o 

-  -  - - 

1 

:  -  o 

o 

i 

1  -  - 

1  -  -  -  - 

0 

. ..  . .  '  I 

!  o 

O' 

o 

. ,o 

o 

Velocity 

magnitude, 

m/sec 


14.0 

13.0 

12.0 

11.0 

10.0 

9.0 

8.0 

7.0 

6.0 

5.0 

4.0 

3.0 

2.0 

1.0 

0 


Figure  1 7. — Computational  flow  field  map;  Reynolds  number,  45,000. 
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1 3.  ABSTRACT  (Maximum  200  words) 

Visual  and  quantitative  information  were  obtained  on  heat  transfer  and  flow  in  a  branched-duct  test  section  that  had 
several  significant  features  of  an  internal  cooling  passage  of  a  turbine  blade.  The  objective  of  this  study  was  to  generate 
a  set  of  experimental  data  that  could  be  used  to  validate  computer  codes  for  internal  cooling  systems.  Surface 
heat  transfer  coefficients  and  entrance  flow  conditions  were  measured  at  entrance  Reynolds  numbers  of  45,000, 335,000, 
and  726,000.  The  heat  transfer  data  were  obtained  using  an  Inconel  heater  sheet  attached  to  the  surface  and  coated  with 
liquid  crystals.  Visual  and  quantitative  flow  field  results  using  particle  image  velocimetry  were  also  obtained  for  a  plane 
at  midchannel  height  for  a  Reynolds  number  of  45,000.  The  flow  was  seeded  with  polystyrene  particles  and  illuminated 
by  a  laser  light  sheet.  Computational  results  were  determined  for  the  same  configurations  and  at  matching  Reynolds 
numbers;  these  surface  heat  transfer  coefficients  and  flow  velocities  were  computed  with  a  commercially  available  code. 
The  experimental  and  computational  results  were  compared.  Although  some  general  trends  did  agree,  there  were 
inconsistencies  in  the  temperature  patterns  as  well  as  in  the  numerical  results.  These  inconsistencies  strongly  suggest  the 
need  for  further  computational  studies  on  complicated  geometries  such  as  the  one  studied. 
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